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On Certain Possible Abbreviations in the Computation 

of the Long-Period Inequalities of the Moon's 

Motion due to the Direct Action of 

the Planets. 

By G. W. Hill. 



Hansen has characterized the calculation of the coefficients of these 
inequalities as extremely difficult. However, it seems to me that, if the shortest 
methods are followed, there is no ground for such an assertion. The work may 
be divided into two portions, independent of each other. In one the object 
is to develop, in periodic series, certain functions of the moon's coordinates, 
which in number do not exceed five. This portion is the same whatever planet 
may be considered to act, and hence may be done once for all. In the other 
portion we seek the coefficients of certain terms in the periodic development of 
certain functions, five also in number, which involve the coordinates of the earth 
and planet only. And this part of the work is very similar to that in which the 
perturbations of the earth by the planet in question are the things sought. And 
as the multiples of the mean motions of these two bodies, which enter into the 
expression of the argument of the inequalities under consideration, are neces- 
sarily quite large, approximative values of the coefficients may be obtained by 
semi-convergent series similar to the well-known theorem of Stirling. This 
matter was first elaborated by Oauchy,* but, in the method as left by him, we 
are directed to compute special values of the successive derivatives of the 
functions to be developed. Now it unfortunately happens that these functions 
are enormously complicated by successive differentiation, so that it is almost 
impossible to write at length their second derivatives. Manifestly then it would 

* Memoire sur les approximations des fonctions de tr&s-grands nombres ; and Rapport sur un Memoir e 
de M. Le Verrier, qui a pour objet la determination' d'une grande inegalite du moyen mouvement de la 
planite Pallas : Comptes Rendus de V ' Academic des Sciences de Paris, Tom. XX, pp. '691-726, 767-786, 
825-847. 
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be a great saving of labor to substitute for the computation of special values of 
these derivatives a computation of a certain number of special values of the 
original function, distributed in such a way that the maximum advantage may 
be obtained. This modification has given rise to an elegant piece of analysis. 
It will be noticed that, in this method, it is necessary to substitute in the 
formulae, from the outset, the numerical values of the elements of the orbits of 
the earth and planet. There seems to be no objection to this on the practical 
side, as, for the computation of the inequalities sought, no partial derivatives of 
B, with respect to these elements, are required. 

I. 

If the masses of the moon, earth and the planet considered are denoted 
severally by m , M and m", and the geocentric rectangular coordinates of the 
moon by a;, y, and z, the similar coordinates of the sun by x', y' and z', and the 
heliocentric coordinates of the planet by x", y" and z", the perturbative function, 
for the direct action of the planet on the moon, is 

p _ »r__ l _{^+^+(y''+y')y+i^+l)n 

But, by a slight substitution in and modification of this expression, we take 
account of the lunar perturbations of the solar coordinates. Let X, Y and Z 
denote the coordinates of the sun referred to the centre of gravity of the earth 
and moon, we shall then have 

x'-X 1 m x i/-Y + ~-^—v r'-Z I m z 
X ~^^ M+m X ' y~ X ^ M + m y > Z -^ + M+m^ 

And A may denote the distance of the planet from the centre of gravity of the 

earth and moon, so that 

A 2 = (x» + Xf+ (*/'+ F) 2 + (z" + Zf, 

also r the radius vector of the moon, so that 

r 2 = x 2 + ?/ + 2 2 ; 

moreover, for brevity, put 

P={x" + X)x + (y"+Y)y + (z" + Z)z. 
Then R takes the form 



B = m" 



T Jf + TO 



M n , M 



l A '-*J^ p + flfW*]' t A! + 2 W^, p + (* T -J"*]'. 
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But it is evident that this expression, differentiated with respect to the 
variables x, y and z, will not furnish differential coefficients identical in value 
with those the expression gives before the transformation, as a/, y 1 and z 1 have 
now been made to involve x, y and z. But a little consideration shows the 
modification which will remedy this. It is plain we ought to multiply the first 

term by *f "* ' and, multiplying the last term by — — . substitute unity 

3 1 
for the numerator and reduce the exponent of the denominator from — to -^ • 

Thus the proper form of B is 



R = ni" 



M r„ ~ M „ , M 



p. M' ,-la 



M+m ' (M+mf 
M-j-m 



+ 



m r * .> , ~ m t^ , «i 2 «~|i 



When this expression is expanded in a series proceeding according to ascending 

powers of the lunar coordinates, and the terms independent of the latter omitted, 

we get 

., ( 4.3 P 2 2 2.1 r* 



2.4 J s 1 2.4 J 3 

M 3 — m 2 r 6.5.4 P 3 __ _3_ 4.3.2 Prh 
+ {M+mf L2X6 l/ r T ' 274^6 1 5 ~J 

i»f 3 +m 3 r8.7.6.5 P 4 _ 4 6.5.4.3 PV 4.3 4.3^2 J. / -1 
+ (Jf + m) 3 L2T4T6T8 ~d»~ "" T'2.4.6.8 _ J 7 + O" 2.4.6.8 J 5 J 

+ 

The terms of this series follow a quite evident law, and it is easy to write as 
many as there may be occasion for. But, hitherto, no sensible inequalities 
have been found arising from the terms beyond the first line. This line furnishes 

all the inequalities which are not factored by the small ratio ■— , whose value is 

about -777- • And the following two lines of terms can add to the coefficients of 
400 

a 2 
these only parts which have the very small factor — li - For these reasons we can 

restrict ourselves to the first line of terms, and write very simply 

3 P 2 1 »- 2 



*=*"L-2 ■■■>— - 2 - A- 
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Restoring the equivalent of P, 

n m (La ^ 5 2 J 8 J xi L2 j 5 2 ^J y 

r 3 y +zy 1 1 -I W + XW+ Y) 

I" L 2 J 5 2 iJ + ? ^ 

+ 3 sww,, + 3 (i+zf+n yz |. 

This expression has the advantage of exhibiting the value of J? as a sum of 
terms of which each is the product of two factors, one of which depends solely 
on the coordinates of the moon and the other is independent of them. 

If we denote the factors of sc 2 , y % and z 2 in R severally by A , B and G, 
we shall have the relation A + B + (7=0. 

Hence it is plain that the number of terms can be reduced from six to five. As 
we shall take the ecliptic for the plane of xy, we will have Z=0. "We can 
then write 

* = HT[7-'-?](«-^ + T y + 17 7 ( * +X V - rf ) 

II. 

We will now express the five factors of the terms of R, viz. r 2 — 3s 2 , 
as 2 — y 2 , xy, xz, and yz, as functions of t, the time, when elliptic values are 
attributed to the coordinates, leaving, however, the longitudes of the perigee and 
node indeterminate, so that the latter may have their motions proportional to t. 
Using Delaunay's notation, and, in addition, putting v for the true anomaly, 
we have x = r cos (v + g) cos A — (1 — 2y 2 ) r sin (v + g) sin A , 

y = r cos (v + g) sin A + (1 — 2y 2 ) r sin (v + g) cos h, 
z = 2y s/l — y 2 r sin (v + g) ; 
or, in a slightly different form, 

x = (1 — y 2 )r cos (v + g + A) + y % r cos (w + g — h), 
y = (1 — y 2 )r sin (t? + g + /») — y 2 r sin (t> + <7 — A) , 
2 = 2yVl — y 2 r sin (v + g)- 
From these equations we derive 
z 2 = 2y 2 (l — y»)r»[l— cos 2(t> + #)], 
r 2 — 3z 2 = [1 — 67 s + 6y 4 ] r 2 + 6y 2 ( 1 — y 2 ) r 2 cos 2 (t; + y) , 
x* — if = ( 1 — y 2 ) V cos 2 (t) + g + A) + 7V 2 cos 2 (» + </ — 7t) + 2y 2 ( 1 — y 2 ) r 2 cos 2A , 
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2xy = (1 — y 2 )¥sin2(» + g + h) — yVsin2(v + g — h) + 2y 2 (l — y 2 )r 2 sin 2A, 
as = y(l — y*)Vsin (2» + 2g + A) + y 3 (l — y 2 )Vsin (2t> + 2gr — A) 

— y (1 — 2y 2 )(l — y 2 )V sin h, 
yz ——y(l — y2)V cos (2t> + 2^ + A) + y 3 (l — y 2 )V cos (2t> + 2# — A) 

+ y(l — 2y 2 )(l — y 2 )VcosA. 
It is then plain that the development of these five factors depends on that 
of the quantities r 2 , r 2 cos 2v and r 2 sin 2». Denoting the eccentric anomaly by 
u, we have ^ = (1 _ e cos tt)2> 

r 2 3 / 1 \ 

-j- cos 2« = -^-e 2 — 2ecos« + ( 1 r-e 2 ) cos 2w, 

— rsin 2» = vl — e 2 (sin 2w — 2e sin «) . 
The constant terms of these functions, in their development in periodic series 
involving multiples of the mean anomaly, are the same as the constant terms of 

the right members of the last equations after they have been multiplied by 

3 5 
1 — e cos u. That is, these terms are severally 1 + -^-e 2 , ye 2 and 0. To obtain 

the remaining coefficients, we put s = s uV ~ 1 , and z = e iv ~\ and recall the theorem ' 
that the coefficient of z i , in the development of any function S according to 
powers of a , is the same as that of s f in the development of 

i ds ' 

according to powers of s. Moreover, adopting Hansen's notation for the 

Besselian function, we put ^'~*'=. 2 4 . /i'V, 

so that, for positive values of i, we have 

dK ~1.2...il l.(» + l)^1.2(» + l)(» + 2) •"J' 
and, for negative values, J^-^ = J^{. 

These functions satisfy the following equation, 

Whence 

.i ,*~ "■' ~ — 

; 



r(f-l) — _ Td) 7(*+D 

"k — , "k "k > 



and, by writing i — 1 for i in the first of these and i + 1 for i in the second, 

Jx («-» = L=l /x («-i)_ / «) f 

J(i + 2) _ «_±1 Jit + 1) _ JW 
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Consequently jmi.^+^I^.^w-ii.^!)^!!], 
The coefficient of z i in the expansion of — being equal to that of s* in 

-K'-K'+tH-t)'* 1 " 8 . 

iS « f r(i-l) x« + l) e / 7 (*-2) r (i + 2) M 

A \" te - " «• <T V " '« " ie 1 \> 

* u 2 2" A \ -jr -j /_! 

which, by means of the relations between the /functions just given, reduces to 

__ JL r (f) 

* 2 

Hence we have ^ „ ^^ . 

_ = 1 + t 6»-2 ii/ .- 7 J|cOB«. 

i=l 

This result may also be obtained from the equation 

r_5l = 2 2. 

<# 2 r 

In like manner we get 

^oos 2t , = 4^+g.4[(l-l«')(/|-'W|«)- e (/|-"-4 + »)]co S ,7, 

i = l 

4- sin 2« = vl=? V .4 W- 2) + ^ +2) - e (/r "+ X i+1 0l sm «. 
a r /• iu-j- -j- NT" T/- J 

i = l 

Consequently, if we put 

+ eB i„»|.7| +,i -(sm'|— i^j-r- 1 ], 

where sin <?> = e, and we agree that 

we shall have, a denoting any arbitrary angle, 

r 2 cos (a + 2v) = a 2 > . if (i) cos (a ~\-il), 

i = — oo 

^sin (a + 2v) = a 2 J^ . H w sm (a + il) . 

l = — 00 

We can now write the expansions of the five factors of the terms of R 
which depend solely on the moon's coordinates : 
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+ y/(l -r 2 )2.£- (i) cos (2g + t7), 
3 - a; '~~ ^ - - (1 — yj 2 . H^cos (2A + 2^ + «) 



4 



1 tW. 



- 3y 2 (1 - y») 2 . 4- «C cos (2A + #) 
+ |~y 4 2.5" (<) cos (— 2A + 2gr+ «), 
4^= l(l-y 2 ) 2 2.fi" (i) sin(2^ + 2^ + *7) 



2 a 2 4 

*" T 
— -J-y 4 2.# (1) sin (— 2A + 2# + if), 



— 3y 2 (l — y 2 ) 2.-4" Jj? sin (2ft + tf ) 



4 

T-5= -|y(l-y 2 ) f 2.^>sin(A + 2^ + *7) 

+ 3y ( 1 - 2y 2 )( 1 - y 2 )* 2 . 4" ^ sin (h + il) 
+ |-y 3 (l - y 2 )*2. JT w sin (- h + 2g + »7) , 

¥? = -4y( i -^ l2 - j3 ' li)cos ^ +25r+ * 7) 

— 3/ ( 1 — 2 r 2 )( 1 — ff 2 . 4 Ji fl cos {h + il) 

+ -|Y(l-y»)*2..ff«C0B(— A + 2flr + »7). 

The summation must be extended to all integral values, positive and negative, 
zero included, for i. When i=0 we must suppose that -^ J '_««_ takes the value 

It will be perceived that the three first terms of R furnish inequalities 
whose arguments do not involve the longitude of the moon's node or involve it 
in an even multiple. The two remaining terms furnish inequalities having an 
odd multiple of this longitude in their arguments. And it is evident that these 
statements remain true even when the solar perturbations of the lunar coordi- 
nates are taken into consideration. Hence, in deriving any particular inequality, 
we never have to consider more than three out of the five terms of R. When 
we propose to neglect the solar perturbations, it can be seen at a glance what 
terms of the expressions above ought to be retained. Thus, in the case of 

Vol. VI. 
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Hansen's inequality of 273 years, the argument involving only I without either 
h or g , it is plain that the first term of — j-^— can alone furnish it ; and, conse- 
quently, we may put, very simply, 

B = _ W V(1 - 6f+ 6yV| [1. - 3 ■£] cos I. 

And the whole difficulty is reduced to finding, in the development of 

1 i* 

the terms A (e) cos (181"— 16Z') + A M sin (18Z"— 16Z') . 

III. 

We pass now to the consideration of the development, in periodic series, of 
the factors of the terms of B which depend on the coordinates of the earth and 
planet. Let it be required to discover the. coefficient C i% v of zV in the develop- 
ment of any periodic function of the eccentric anomalies u and v! of two planets, 

in the case where i is quite large. We shall suppose that the function has -^ 
for a factor. It is known that 

_L = N m j-j _ 2a cos (m - §) + a 2 ]"" [1 - 2fr cos (u+Q) + W]~ n , 
where N, a, ft and Q are functions of u! or V only, and VC and ft are always less 
than unity. Substituting the imaginary exponential s=s uV ~ 1 , and, to abbre- 
viate, putting & = ar 1 e«> / - 1 , k 1 =Jr 1 £- QV ~ 1 , 
this equation becomes 

1 r- J , we can then suppose that the function 

to be developed is 

(i -!)">). 

The coefficient of z l in the development of this is equivalent to 

<*= if-'-* ( - •» t 1 - \ (• + 1)] (> - Tpw *■ 

Let us put /(<) = ,4<-4) [i _ JL (, + 1)] ^.) . 

Since the absolute term of a series of integral powers of a variable is not 
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changed by substituting for the latter a constant multiple of it, in the expression 
for G t we can write fos for s. Thus 

The difficulty here that the factor (1 — s)~ n becomes infinite at the limits of the 
definite integral, is only apparent. For the multiple of s instead of 7cs may be 
ps, in which the modulus of p is less than that of h by a very small quantity. In 
this case we get a tangible result, which is seen to have, as its limit, when p is 
made to approach Tc indefinitely, the value which will be presently given. 

We now assume that it is possible to expand f(Jcs) in an infinite series 
proceeding according to positive integral powers of u* Let 

f(ks) = C -|- CjM + C 2 W 2 + . . . =%.CjU J . 

Then (7,= |^2. J^l-^- 1 (1 — e^-^c^du. 

The definite integral J_ j^-tuv-in _ ^v-iy—Ou 

is a function of n and i: with Oauchy we will denote it by [n} t . Then by taking 
the derivative of the quantity, under the integral sign, j times with respect to i, 

we S et i_ .J^" 6 -iuv-i (1 _ E uv-i)-n u j du _ ^zzrjyDi . [n] t . 

Whence we have the symbolic expression for G t , 

C^k-ifike-^.lnl. 
But We have «».= 1 + A, <-».= jA., 

A here denoting the characteristic of finite differences with respect to the 
variable i , and not the distance between the two planets. Let 

V= 1 I j > then e~ d '=1 — v- 

Making these substitutions, we have 

C7 4 = *-'/(&-&v)-[«]i- 
By successive integrations by parts, making the integration always bear on 

the first factor, we find the value of the definite integral, 

^/"'-"-d- *■-)-*.= w,= " ( " + \\i ; ^t- 1 ' ■ 

*This is the assumption which leads to the semi-convergent series representing the value of Q. 
Its allowableness is shown by the fact of the relative smallness of the definite integral which ought to 
be added to complete the truncated series, when i is tolerably large and the number of terms taken 
into account is not too great. As Cauchy has treated this point at length, in his memoir first mentioned 
above, I have thought it unnecessary to say more about it here. 
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When the function f(h — fay) is developed in ascending powers of v » the general 
term of G t will be proportional to 

A } 

And, developing the last expression for G t , and employing accents, attached to /, 
to denote differentiation of the form of /, we have 

G i = h-^f{h)[n^~hf{h)in~l^+~¥f'{h)\n-1i\ i 
This may also be written 



C= *-•[»], | /(*) -fW-ij^i 



+ -L rW .J._>- 1 X»-2) 



l^'' v '' (,' + n — l)(» + n — 2) 

1 -fi(]c).tf^ (»-l)(n-2j(»-S) 



1.2.3- 7 v y (j + »_i)(i + n _2)(» + n — 3) 

+ 1 

We may employ the T function to 'express [n]*, and then 



/»/■(* +i) 

In practice, n will have some one of the following series of values, 

1 3 5 7 

T' T' T' T* etc -' 

and it is well known that 

r (!) = </*, r (4) = 1r^> r (y) = T V7t - etc - 

When « is a tolerably large integer, we may use the semi-convergent series 
log r (t + n) = y log (27t) +(t + n — -y) log (» + n — 1) 

logr(»+l) = ilog(27t) + (i + ^-)log* 

+ iK ( * + 1.2 i 3.4 t^ + 5i? ••• J ' 

where Mis the modulus of common logarithms, and B lt B 3 , etc., are the numbers 
of Bernoulli; Thence is derived 
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^ r ^ ) ={i + \yo g i ±^ + {n-l)\og{i + n-l) 

_ uS i ^ f 1 1 -] .B 8 r 1 1 -i 

M \ n + 1.2li i + n — J 3AL? (i + n — 1) 8 J 

+ Afl _J I- 1 

^5.6 L* (i + n—l^J ••■/ 

= (i + \)\og i + n .- 1 +(n-l)\og(i + n-l) 



—M{n — 1) | 



1 + 



1 1 1 i» + »(t + n— !) + (» + n — l) 2 



+ 



12 »(< + n — 1) 360 j»(i + n_l)8 



1260 »*(i + B _ 1)' 



,6 



The first term of the last expression for G { affords a first approximation to 

its value, correct, so to speak, to quantities of the order of — • Then 

£, = &-<[>],/(&). 
In like manner, the two terms at the beginning afford an approximation correct 

to quantities of the order of -^-- Here we can effect a remarkable reduction; 

for, on comparing the two terms in question with the two first terms of Taylor's 
theorem, we see that, to the same degree of approximation, we may write 

No more labor is involved in employing this expression than in the preceding. 

IV. 

In this condition Cauchy leaves the subject, but we may go a step farther. 
In the cases which come up in practice f(k) is always such a function that 
successive differentiation immensely complicates it ; so that it is scarcely possible 
to go beyond f"(k). Hence a great deal of labor is saved, if, instead of 
attempting to calculate /' (k) , f" (k) , etc., we substitute the calculation of f(k) 
for several values of the argument k. It is easy to perceive that, in general, all 
the derivatives /' (k) ,/" (k) , etc., may be eliminated from the expression for G t . 
For, cutting the series off at the term which contains f m (k) as a factor, we may 
suppose that, to the same degree of approximation, 

G i = k~ i [n]« { x f{k — ky ) + x x f{k — ky^+... +x p f(k — ky p ) } , 
where x , £c 1? . . . x p and y , y lt . . . y p are unknowns to be suitably determined. 

By developing this expression for G t in powers of k and comparing it with 
the previous expression, we get the following system of simultaneous equations 
for determining the unknowns x 0) a^, . . . x p , y , y u . . . y p : 
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*o "f »i + %2 H- . • . +x p = 1 , 

a<>2A>+»a2/i+Z22/2+- ••+«!>%> = j ? ft _i ' 

.«WlJ + «iJ« + a?i»8 + ...+a»!4 = (i +*Z 1 ix" + n-2) ' 

~ ? ,2?+i+~ w 2j>+ij_~ w 2p+i_i_ 4-3. ,/ 2 p + 1 = (»— 1).. -(n — 2p — 1) 

z<>2/o i-«i2/i -1-a»y» + '-' + J »& (» + n — 1)... (i + n — 2p — 1) 

For the sake of brevity we will denote the right-hand members of these 
equations as « , a lt a % . . . a ip+1 . The solution of these equations is very 
elegant. According to the theorem of Bezout, the degree of the final equation, 
obtained by elimination, would be = (2p+ 2)! But, as we shall see, the solu- 
tion depends on that of an equation of the (p + l) th degree, whose roots are the 
values of the several unknowns y , y x . . . y p ; and there is practically but one 
solution. 

Let us suppose that the values of the y's, in any particular solution, are the 
roots of the equation 

y p+1 + siy p + s 2 y p - 1 + ...+s p+1 = o, 
so that *i= — (2/0 + 2/1 + 2/2 + • • • +%,). 

% = 2/o2/i + 2/n2/a + 2/i2/2 + • • • , 



s p+1 = (— 1) p+1 2/o2/i2/2 • • • 2/j>- 
Hence, y q denoting any one of the j/'s, we must have 

2/? +1 + «i2/r + ny?- 1 + • • • + s p+i = ° • 

Now, in the group of equations to be solved, multiply the equation, whose 
second member is a p+u by 1, the one, whose second member is a P , by s lt and 
so on until the first equation is multiplied by s p+1 . Then, by adding all the 
equations thus obtained, the first member of the resulting equation vanishes, and 

we have a P + i + a p Si + a P -i s 2 + • • • + «o- s i>+i — °- 

By cutting off the first equation and adding to the group the equation whose 
second member is a p+i , and writing a" c for x y , Xi for x^y-^, and so on, we obtain 
a group which differs from the former only in the second members. Hence we 
have, from this group, the equation 

<ip+t + a P +iSi + <t' P s 2+ • • • +«i*p+i=0. 
And, in a similar manner, 

#p.f3 + #p+2 s i + a p+i s 2 + • • • ~r %#p+i == 0> 



«2j>+i + «2p*i + < l 2p-i s 2+ • • • -\-a p s P + i — 0. 
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These p + 1 linear equations suffice to determine the values of s l: s 2 . . .s P+1 , 
the coefficients of the equation of the {p-\- l) th degree, which has, as its roots, 
the values of the unknowns y. These values, being obtained and substituted in 
the first {p + 1) equations of the original group, we have a group of (p + 1) 
linear equations for determining the (p + l) unknowns x , a- x . . . x p . It is plain 
that all possible solutions of the group of equations are obtained by permuting 
be-tween themselves the roots of the equation which gives the values of the y's ; 
and, as thus, to each root, corresponds its special value of the x's, and the order 
in which the several terms of (7 4 stand, is of no import, it is clear that, practically 
at least, but one solution exists. 

In practice, p never need exceed 2. For p = , the solution has already 
been given. For p = 1 , we have 

(»-lXn-2) n-1 

(n-l)(n— 2)(n-3) {n -l)(n-2) n-1 

(i + n —l){i + n—2){i + n— 3y {i + n — l){i + n— 2) 1 ^ i+n— 1 2 
The solution of these gives 

_ n n — 2 _ (n—l)(n — 2) 

*i— 'i + n —3' S *~ (i + n — 2)(i + n — 3)' 
Thus the equation which contains the values of the y's is 

, n—2 (n— l)(w— 2) _ 

y 2 i + n ^3^ + {i + n — 2){i + n — S)~°' 

Whence the two values of y are 



„ , /(2— w)(t— 1) 
n — 2 ±v/ > — , - — 5-^ 

V— i + n — 3 ; 

and the corresponding values of a are 



_ 1 r i — n + 1 / i + n — 2 ~| 
• a3_ 2 L 1 ± t + n— 1 V (2 — n)(t— 1)J ' 
In many cases these values will be imaginary, which, however, does not hinder 
their use, as h is imaginary. 
For p = 2 , we have 

(n—l)(n— 2){n— 3) (n— l){n — 2) n— .1 , 

(i + w _l)(i + n _2)(t + w — sy {i + n— l)(» + n — 2)* 1-| "i + n— l Sj + * 3 

(w _2)(n-3)(n-4) + (n-^-S) n-2 _^ 



(i + n — 2)(i + n— 3)(i + n — 4)^ {i + n — 2)(i + n — 3) x ' t + n— 2 
(„_3Xn-4X»-5) ^ (»_3Xn-4) n-3 



({ + »_3)(«- + »— 4X* + n — 5) ' {i + n — 3,){i + n — 4) x ' t + n— 3 
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The solution of these equations gives 

. - ■> n ~ S »-o (n-2)(n-3) (n-l)(n-2)(n-3) 



t+ n _5 * ~ (,-+„_ 4)(i+ n _5)' "' - (,-+»_3)(i+ n _4)(<+»_5) 
The equation, which has, for its roots, the values of the y's, is 
?J 3 , "~3 ... (n-2)(n-3) (n-l)(»-2X»-3) _ 

* °; + „_5^ T °(i + n — 4)(« + n — 5) y •(« + »» — 3)(i + n— 4)(i + n — 5) 
By comparing this with the equation for the case where p = 1 , we readily see 
what the equation would be for higher values of p. 

As an example, suppose it were required to find the coefficient of z 18 in the 
expansion of [l — 2a cos (u— Q)+ u 2 ]~ *. 

Here the form of /(«) is 

In the first place let two terms in the final expression for C t be regarded as 

3 
sufficient, that is, put p = 1 . Then i = 18 , « = -»-, and the two values of y are 

33 
and the corresponding value of a; is 



■ = T('*^Vf,)- 



Thus the expression for C< is 

C 18 = k- 1B [~~\ {1.17865/(0.9880647%) — 0.17865/(1.0725413%)}. 

L I Ji 8 

The error of this is of the order of -^ , while, in case p = , which gives the 

which Cauchy employed, the error is of the order of -^-- 

In case we make p = 2, and thus have three terms in the formula for G t , 

the roots of the cubic 3 , 9_ 2 . 9 9 

y ~*~ "29 y "*" 31.29 2/i " 33.31.29 ~ 
must be found. They are 

?/ = + 0.00804343, y, = — 0.04617994, y % — — 0.27220828. 

The linear equations for determining the x's are 

x + a?i H- «2 = 1 , 

0.0804343 a: — 0.4617994^— 2.722083x 2 = 0.2702703, 
0.0064697a: + 0.2132586a! + 7.409736^ = — 0.0772201. 
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The solution of which gives 

£C = + 1.3426685, a^ = — 0.3408857, x % — — 0.0017828. 
Thus, in this case, we should have 

C 19 z=h- U \~~\ S 1.3426685/(0.9919566%) — 0.3408857/(1.046 179947c) 

— 0.0017828/(1.2722083%)}. 
The error of this formula is only of the order of -^ • 

In further illustration of this method, let us find the value M 8) of the coeffi- 
cient of cos 180 in the periodic development of 

(1 — 2a cos + a 2 ) -1 , 
where a = 0.723332 the ratio of the mean distances of Venus and the earth from 
the sun. Here the form of f{s) is simply 

/w=(>-fr 

Let us take the formula where p = 1 . We have 
*f=2tf 18 =2[4] i8a «j 1.17865(1^ 

The value of f-^-1 will be found in the table at the end of this memoir. And, 

on the substitution of the numerical values, we get M 8) = 0.090880. Delaunay, 
in his memoir,* has 0.090876. 

In the case where the function to be developed contains the anomalies of two 
planets, after the value of C t has been obtained corresponding to j points evenly 
distributed on the circumference with reference to the variable V or the variable 
u', the value of C t v results by employing the method of mechanical quadratures : 
the formula in the first case being 

C i<i ''=Z—r'2<.C i Z~ , 

and, in the second, n \ -. n r 1 ,_ v 

i,v ~J * *~d ' 

In the annexed table are given the common logarithms of the function [rijt, 
9 
for n as far as n = -=-, and for i, as far as i = 30. As they have been computed 

with the ten-figure logarithms of Vega's Thesaurus Logariihrnorum, it is to be 
presumed that they are correct, in nearly every case, to half a unit in the last 
place. 

*Connaissance des Temps, 1862. 
Vol. VI. 
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Table op the Values op Log [n] ( . 



i. 


» = i- 


n = f. 


« = #• 


» = *• 


n = h 


1 


9.6989700 


0.1760913 


0.3979400 


0.5440680 


0.6532125 


2 


9.5740313 


0.2730013 


0.6409781 


0.8962506 


1.0925452 


3 


9.4948500 


0.3399481 


0.8170693 


1.1594920 


1.4283373 


4 


9.4368581 


0.3911006 


0.9553720 


1.3703454 


1.7013386 


5 


9.3911006 


0.4324933 


1.0693154 


1.5464366 


1.9317875 


6 


9.3533120 


0.4672554 


1.1662254 


1.6977043 


2.1313599 


7 


9.3211273 


0.4972186 


1.2505463 


1.8303299 


2.3074511 


8 


9.2930986 


0.5235475 


1.3251799 


1.9484292 


2.4650590 


9 


9.2682750 


0.5470286 


1.3921267 


2.0548845 


2.6077265 


10 


9.2459986 


0.5682179 


1.4528245 


2.1517945 


2.7380602 


11 


9.2257953 


0.5875231 


1.5083418 


2.2407356 


2.8580356 


12 


9.2073118 


0.6052519 


1.5594944 


2.3229224 


2.9691860 


13 


9.1902785 


0.6216423 


1.6069190 


2.3993107 


3.0727266 


14 


9.1744842 


0.6368822 


1.6511227 


2.4706666 


3.1696366 


15 


9.1597610 


0.6511227 


1.6925154 


2.5376134 


3.2607171 


16 


9.1459727 


0.6644866 


1.7314334 


2.6006651 


3.3466317 


17 


9.1330077 


0.6770758 


1.7681562 


2.6602508 


3.4279367 


18 


9.1207733 


0.6889750 


1.8029183 


2.7167322 


3.5051026 


19 


9.1091914 


0.7002560 


1.8359186 


2.7704171 


3.5785315 


20 


9.0981960 


0.7109799 


1.8673271 


2.8215696 


3.6485694 


21 


9.0877306 


0.7211990 


1.8972903 


2.8704181 


3.7155162 


22 


9.0777464 


0.7309589 


1.9259355 


2.9171615 


3.7796337 


23 


9.0682010 


0.7402989 


1.9533737 


2.9619739 


3.8411517 


24 


9.0590577 


0.7492537 


1.9797027 


3.0050085 


3.9002732 


25 


9.0502837 


0.7578539 


2.0050085 


3.0464012 


3.9571780 


26 


9.0418506 


0.7661264 


2.0293679 


3.0862727 


4.0120267 


27 


9.0337327 


0.7740954 


2.0528490 


3.1247310 


4.0649628 


28 


9.0259073 


0.7817822 


2.0755129 


3.1618728 


4.1161153 


29 


9.0183542 


0.7892062 


2.0974148 


3.1977853 


4.1656007 


30 


9.0110560 


0.7963858 


2.1186051 


3.2325484 


4.2135252 



